Decoherence and Full Counting Statistics in a Mach-Zehnder Interferometer 
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We investigate the Full Counting Statistics of an electrical Mach-Zehnder interferometer pen- 
etrated by an Aharonov-Bohm flux, and in the presence of a classical fluctuating potential. Of 
interest is the suppression of the Aharonov-Bohm oscillations in the distribution function of the 
transmitted charge. For a Gaussian fluctuating field we calculate the first three cumulants. The 
fluctuating potential causes a modulation of the conductance leading in the third cumulant to a 
term cubic in voltage and to a contribution correlating modulation of current and noise. In the high 
voltage regime we present an approximation of the generating function. 
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I. INTRODUCTION 

For two decades theoretical and experimental inves- 
tigations of the noise properties of small electrical con- 
ductors have been an important frontier of mesoscopic 
physics 1 . Non-equilibrium noise, in particular shot 
no j se 2 J 3 J 4 J £ can p rov ide information on physical proper- 
ties which are not accessible from conductance measure- 
ments. Recently there has been considerable interest in 
characterizing transport in mesoscopic systems not only 
via conductance and noise but also with all higher or- 
der current correlations. This is achieved by deriving 
a generating functional from which the quantities of in- 
terest can be obtained simply by taking derivatives. To 
construct the generating functional, it is useful to imag- 
ine that electrons passing the sample in a given time 
are counted, and the approach is thus known as Full 
Counting Statistics^. While initial work was based on 
the scattering approach to electrical transport, subse- 
quently a number of different methods have been de- 
veloped. These include an approach based on Keldysh 
Green's functions&2ii2iii, the non-linear sigma modeliS, 
an approach based on a cascade of Boltzmann-Langevin 
equations^ and a formulation in terms of a stochastic 
path integral 14 ' 15 ! 1 "' 1 ^ 18 ! 19 ! 2 "! 21 ! 22 . These later methods 
are principally useful if the quasi-classical part of the 
transport coefficients is dominant. A first experiment 23 
has given additional impetus for this research. It has 
stimulated work which treats the effect of the measure- 
ment circuit on the counting statistics^ and led to pro- 
posals for novel detection schemesS^SiSi. 

In this work we are interested in the Full Counting 
Statistics (FCS) of transport problems in which quan- 
tum effects are important. For instance, in single channel 
mesoscopic conductors quantum interference is essential. 
In such a situation, transport becomes classical only due 
to dephasing. It is thus interesting to examine the effect 
of inelastic processes on the counting statistics of a single 
channel problem and to investigate the transition from 
quantum mechanical transport to fully classical trans- 
port. 

The central aim of our work is thus the derivation of 
the generating functional in the presence of inelastic scat- 




FIG. 1: The Mach-Zehnder interferometer: a four-terminal 
AB-ring, threaded by a magnetic flux. The interferometer 
arms are coupled to classical fluctuating potentials modeling 
the effect of dephasing. The influence of dephasing on current 
fluctuations in lead 3 is investigated. 



tcring. We investigate an interferometer penetrated by 
an Aharonov-Bohm flux and subject to a random classi- 
cal fluctuating potential. For simplicity we assume that 
the interferometer is of the Mach-Zehnder (MZI) type 
(Fig. nj, an interferometer without backscattering. In 
such a system there are no closed orbits and thus only a 
finite number of possible trajectories. 

Despite these simplifying assumptions the derivation 
of the generating functional is a non-trivial problem. To 
investigate the effect of the fluctuating potential, a time- 
dependent treatment of transport is needed together with 
a technique of statistical averaging. In fact, in the present 
work we are able to give the generating functional only 
in the limiting case of the high voltage regime. In the 
general case, we present results for the first three cumu- 
lants. 

The model we investigate is closely related to earlier 
work by Seelig et al. on the effect of dephasing on the 
conductance^ and four-terminal resistances in ballistic 
interferometers^ 9 ". In few channel structures dephasing 
might be sample specific 29 . Experimental investigations 
on dephasing in ballistic interferometers is provided by 
Hansen et ali^S and Kobayashi et ali^i. Recently an elec- 
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tronic MZI has been experimentally realized by Ji et ali2£, 
using edge states in the quantum hall regime. This exper- 
iment also raised interesting questions on the shot noise 
of the interferometer and dephasing. The effect of de- 
phasing on shot noise in such an interferometer has been 
discussed by Marquardt and Bruder based on classical^ 3 " 
and quantum 34 fluctuating field models. 

Still another motivation for the work presented here 
derives from recent proposals to implement orbital 
entanglemenfeS based on electron-hole generation^ in 
two-particle intensity-interferometers 3 ^. In these propo- 
sals quantum interference effects play an important role 
in noise and higher order correlations 3 ^. 



II. THE MACH-ZEHNDER INTERFEROMETER 



The conceptually simplest system displaying electron 
interference is the electronic Mach-Zehnder interferome- 
ter. The setup is shown in Fig.^ it consists of two beam 
splitters connected by two interferometer arms enclos- 
ing a magnetic flux. Each beam splitter is characterized 
by a transmission amplitude \pT for an electron going 
straight through, and a reflection amplitude iVTZ, with 
1Z + T = 1. The simplicity of the interferometer origi- 
nates from the exclusion of backscattering which discards 
closed orbits and separates completely the processes of 
left and right moving particles. 

Interference occurs, because the different vector poten- 
tials in the two arms lead to a phase difference between 
the paths A and B. This difference creates a charac- 
teristic flux-periodicity in the interference pattern, the 
Aharonov-Bohm effect. 

In our work, dephasing is introduced with the help of 
classical fluctuating potentials, J7a(£) and Us(t). The 
correlation function of these fluctuating potentials can 
be self-consistently determined treating electron-electron 
interactions in a random phase approximation^. Here 
we take this correlation function to be given. The poten- 
tials cause the particles to pick up an additional time- 
dependent phase <fi(t) — ft +T dt'Ui(t'), i = A or B (here 
and in the following we set e = H = 1). We assume, that 
the traversal time r of the particles is not affected by the 
potentials. The transmission amplitudes of particles en- 
tering the interferometer at time t in lead (3 and leaving 
it at time t' in lead a are collected in the time-dependent 
scattering matrbsS S a p(t,t'). For a one-channel MZI 
this scattering matrix is four dimensional and has spe- 
cial properties: as backscattering is excluded, the scat- 
tering matrix has non- vanishing elements only in the off- 
diagonal 2x2 blocks. Because the motion is purely bal- 
listic, the scattering matrix depends effectively on one 



single time S(t, t') = 5{t' - t - r)S(t), where 

s 31 {t) = iVrczV* 4 -^ + 1) (I) 

S 32 (t) = Te^ +Mt) -H 
S 4 i(t) = -fte 4 * +Mt) +T 
S i2 {t) = n/^ffe^+^W + 1). 

Here we assumed equal arm lengths and introduced the 
phase difference tp(t) = <Pa{£) — and the magnetic 

flux $ in units of the flux quantum. The remaining four 
elements are found from the magnetic field symmetry of 
the scattering matrix S a p(&) — Sp a (— $). 

Without the fluctuating potentials the system is co- 
herent: tp(t) = 0. The transmission probability becomes 
T31 = 1 531 1 2 = TR!T{\ + cos<I>) and shows cosine oscilla- 
tions in flux. Such AB-oscillations are visible in conduc- 
tance experiments. 

When we include dephasing by the fluctuating poten- 
tials, the scattering matrix becomes time-dependent and 
the quantities of interest have to be averaged over the sta- 
tistically distributed potential. This average suppresses 
the AB-oscillations. For the case of conductance and 
noise, these effects have been studied in detail 2 ^*^ 3 .. In 
this work we will generalize this discussion to higher or- 
der current correlations (Full Counting Statistics). 

III. FCS OF THE COHERENT SYSTEM 

The Full Counting Statistics^ (FCS) is the distribu- 
tion function P(Q) giving the probability that a certain 
charge Q — f Q ° dtl(t) flows through a system during the 
observation time to- Its Fourier transform is called its 
generating function x(A), the conjugated variable A is 
named the counting field: 

X (A) = 5>(QK AQ (2) 
Q 

P(Q) - ^ f d\ X (X)e~ lXQ . (3) 

The expansion of ln%(A) in A yields all irreducible cumu- 

lants: lnx(A) = Y^k ^Tf - (Q k )- ^ or l ar § e *o the first cu- 
mulant gives the mean current (J) = (Q) /to, the second 
cumulant is the zero frequency noise S = (Q 2 ) /to, the 
third cumulant C 3 = (Q 3 ) /to describes the asymmetry 
(skewness) of the distribution, and so on. The generating 
function provides in principle a complete description of 
the zero- frequency current fluctuations. 

We consider a four-terminal conductor, thus Q and A 
are four dimensional vectors. Since in the MZI current 
conservation relates current fluctuation of all leads, it 
is in most cases sufficient to concentrate on the charge 
distribution in one lead. In the following we will chose 
lead 3 and call Q its charge and A the corresponding 
counting variable. 



In the coherent regime {(pit) = 0), the scattering ma- 
trix is time independent, and the calculation of the gen- 
erating function is straightforward. Following the work 
of Levitov and Lesovik£, one defines an energy-dependent 
4x4 matrix A = f ■ (S* S — 1), where the diagonal ma- 
trix / = f(E) contains the Fermi functions of the four 
terminals. S incorporates the counting fields by phase 
shifts S a p — e ltyXa ~ x ^ S a p (remember that \ a = \5 a3 ). 
The generating function is 



ln X (A) 



2tt 



dEtr\n{\ + A). 



(4) 



Because backscattering is excluded, the matrix A is block 
diagonal, and lnx separates into two independent contri- 
butions. This fact expresses the statistical independence 
of left and right going particles. 

For a two-terminal conductor, the FCS can always be 
interpreted as a classical binomial distribution. For a 
four terminal conductor like the MZI this is in general 
no longer the case. There exists however a special con- 
figuration: at zero temperature a voltage V is applied 
on terminal 1 and zero voltage on the other terminals. 
In this case only two processes exist: transmission from 
terminal 1 into lead 3 with probability T31 or into lead 4 
with probability 1— T31. The number of incident particles 
during the time to is given by N = . We then recover 
the binomial distribution of the two-terminal conductor 



X(X) = (e iX T 31 + (l-T 31 )) N 
P{Q) = (2)^(1 -Tsx) N - Q 



(5) 
(6) 



with the only difference, that the transmission probabil- 
ity T31 = 21ZT(\ + cos<l>) depends on the magnetic flux. 
Consequently the distribution function P(Q) contains all 
harmonics up to cosA^$. The flux periodicity of P(Q) 
is shown in Fig.^a). For symmetric beam splitters and 
certain values of <!>, all particles end up in one lead, due 
to complete constructive interference. The distribution 
is then singular and noiseless: P(Q) = Sq,n for $ = 2nk, 
and P(Q) = S Qfi for $ = (2k + 1)tt with'fc G N . 

At finite temperature this simple binomial picture is 
modified and does not coincide with the two-terminal 
conductor statistics-. We obtain 



X(A) = exp 



2ir 
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arcosh 2 u(A) — v 2 + ivX A z 
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with 



v 



2k B T 



and the complex 



y) 
function 



u(A) = T31 cosh (v + if) + (1 - T31) cosh (y - if) . 

Fig. ffib) shows plots of the logarithm of the distri- 
bution function for finite temperature. At high tem- 
peratures the flux dependence is strongly reduced, as 
is to be expected for a general interferometer. But 
the MZI with arms of equal lengths investigated here 
has the special property that it is described by energy- 
independent transmission probabilities. Thus the con- 
ductance does not suffer from energy averaging and ex- 




FIG. 2: (a) The distribution function at zero temperature 
(logarithmic scale, Vtof2n = 40). Left: 3D-plot of the bino- 
mial distribution, which is periodic in flux. Right: the same 
in contour representation, dark areas correspond to low prob- 
ability, (b) The distribution function at high temperature 
k B T = 10V (for a long observation time t k B T = 100). The 
flux dependence is washed out, and the distribution tends 
towards a Gaussian. 



hibits full Aharonov-Bohm oscillations even at high tem- 
peratures. However, the distribution function becomes 
classical, because the high-temperature fluctuations of 
transmitted charge are dominated by two-particle pro- 
cesses, whose probability is flux independent and equal 
to one. 

At very high temperatures ksT 3> V, we 
find a Gaussian function P(Q) - e -(Q-(Q)) 2 w i t h 
centered around the me 



AQ- 

<Q) = £t 31 v. 



2t k B T> 



charge 



IV. INFLUENCE OF DEPHASING 

We now start to consider the influence of the fluctuat- 
ing potential on the FCS. The scattering matrix becomes 
time-dependent, and the electrons can be scattered be- 
tween different energies. Following the idea of Ivanov and 
Levito\™i, the different energies can be viewed as an infi- 
nite number of independent scattering channels. Switch- 
ing to time-representation S(t,t') = 5(t' — t — r)S(t), 
this leads to a generalization of formula J2J for the gen- 
crating function, where we extend the usual trace (tr) to 
a trace in time-space (Tr) including time integration: 



lnx(A) = Trln(l + A) 



(8) 



All matrices are defined in analogy to sec- 
tion IIIII The matrix A depends on two time- 
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variables A(t', t) = n(t' - t)(S^(t)S(t) - 1); the 
diagonal matrix n contains the Fourier trans- 
forms of the Fermi functions in the four terminals 
n a =n a (t- t') = J |f/ Q (^)e lB (^*'). It can easily 
be checked, that Eq. 10 reduces to Eq. (@J, when A 
depends only on time-differences (elastic scattering). 

To study the influence of dephasing on the gener- 
ating function, it has to be averaged statistically over 
all possible realizations of the fluctuating potential: 

(x(X))u — ^ e ^ rln ( 1+A ^ . For a Gaussian distribution, 

this average is expressed by a functional integral 



(*(A)) = / VUx{\{U})^ h!t ° 0dt!dtlU{t)K{t ^' )U{tl) - 



(9) 

The kernel K is defined via K(t) = J gS^^wJe"^, 
with the spectral function of the potential 
IVH = J dt(U(t)U(0))e iuJt . Only potential differ- 
ences matter: U(t) = [7,4 (i) — £/b(£). 

In general the averaging can not be performed analyt- 
ically. Therefore, in a first step, we expand x(A) in the 
counting field A and calculate the average of single cu- 
mulants, crosschecking known result a 28 i 33 and presenting 
new results for the third cumulant. In a second step, we 
will derive directly an expression for xW m the case of 
a slowly fluctuating potential. 



V. CUMULANTS, EXACT AND LIMITING 
CASES 

The effect of the averaging on the cumulants is deter- 
mined by three quantities: the correlation time trc of 
the fluctuating potential U(t), the traversal time r, and 
a decoherence parameter z. 

trc is typically the i?C-time, with C the coupling ca- 
pacitance, and R the charge relaxation resistance^; it 
provides a cut-off frequency for energy exchanges be- 
tween particles and fluctuating potential. The traver- 
sal time r is the correlation time of the phase tp(t), 
because the electrons integrate the fluctuating poten- 
tial U(t) during time r. For trc <C r the corre- 
lation function of the phase is a triangle in time 42 
(ip(t)ip(0)) = Eu ■ (t - \t\) 9(t - with Sy the spec- 
tral function of the potential at zero frequency. This 
spectral function is given by S;y = AU 2 trc, where AU 
is the mean fluctuation amplitude of the potential. Be- 
cause the relevant quantity in the MZI is the accumu- 
lated phase (p(t), it is convenient to introduce the reduced 
mean fluctuation amplitude AU = y/^SSAU. The spec- 
tral function Sy defines -for a Gaussian fluctuating field 
U(t)- the decoherence parameter z via 2 ^ 3 ^: 



, Syr 
exp( jr-}. 



(10) 



The parameter z indicates the dephasing strength; it 
varies between for complete dephasing and 1 in the 
absence of dephasing. 



FIG. 3: Noise (a), (k B T = 0) and third cumulant (c), 
(ksT — 0.1) as a function of decoherence parameter z and 
flux 4> (for Vt = 1): the AB-oscillations are suppressed for 
decreasing z. (b) and (d) display the dependence on voltage 
and temperature (z = 0.3). The third cumulant changes little 
with temperature, but the noise shows linear thermal noise for 
high temperature. 



We consider the configuration with V applied on ter- 
minal 1, and current fluctuations observed in lead 3. The 
cumulant of fc-th order is given by 



(Q k ) = 



1 dHnjxjX^u 



d\ k 



(11) 



A=0 



The calculation for the first cumulants is sketched in the 
appendix B here we will give only the main results and 
discuss the figures. 

The mean current in lead 3 is linear in voltage, 
and proportional to the transmission probability, aver- 
aged over the potential: (I3) = ^27ZT(1 + z cos &)V. In 
comparison to the coherent system, the AB-oscillations 
in the conductance are suppressed by the decoherence 
parameter z, in agreement with the work of Seelig and 
Biittikes^. 

The noise has been studied in detail by Marquardt and 
Bruder— . The expression for the noise is separated into 
a Nyquist part S nyq , a modulation part Si and a two- 
particle exchange contribution S ox . The ^-independent 
Nyquist noise S nyq describes thermal noised as well as 
energy exchange between particles and the classical po- 
tential U{t) up to the cut-off frequency I/t^c- It is finite 
even for zero temperature and voltage. Si is a conse- 
quence of the variation of the conductance in time, it 
grows with V 2 and is temperature-independent. The ex- 
change term S ex describes shot noise and is a function of 
voltage and temperature, vanishing for V = and linear 
for high voltage. S cx is a consequence of the indistin- 
guishability of identical particles. Thus the total noise is 



1 (T) + E I (V 2 ) + E ex (V,T). 



(12) 
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FIG. 4: Processes contributing to the third cumulant at zero 
temperature. The black points correspond to channels a, /3, 
7 at time t + t, and the crosses to channels p, p, v at time t. 
In Cf (a) three one-particle processes are connected by cor- 
relations of the fluctuating potential (circled in grey/green). 
In (b) the correlation of the modulated current and the shot 
noise are displayed, which contributes to Cf s . C| p (c) is the 
three particle contribution, averaged over excitations. 



For symmetric beam splitters, the fundamental AB- 
period is absent in the noise and the AB-oscillations are 
proportional to cos 2$. They are suppressed for z — > 0, 
compare Fig. Qfa). 

In Fig. |SIb) the noise is shown as a function of tem- 
perature and voltage: it is linear in temperature for high 
fcsTr, simply due to thermal noise. The quadratic de- 
pendence on voltage does not change with temperature. 

In the third cumulant we distinguish three contribu- 
tions: a modulation part Cf, proportional to 7 3 , a sec- 
ond term C 3 S , that connects current modulation and 
noise, and a three-particle contribution C 3 p : 



C 3 = Cf(V 3 



Cf 



is 



(V,T) 



Cl p (V,T) 



(13) 



The processes leading to these three contributions can be 
visualized as shown in Fig. Fig. Ufa) corresponds to 
the modulation contribution Cf. Each circle contains the 
one-particle process creating current in the corresponding 
lead: In the left half of the first circle an electron, entering 
the scatterer at time t in lead p is transmitted into lead a 
at time t+r with the amplitude S ap , it is moving forward 
in time. Correspondingly in the right half circle a hole is 
transmitted from lead a into lead p with amplitude Sj a , 
moving backward in time. We sum over all channels at 
the points X and measure the current at the points •. As 
we investigate the third cumulant of charge transmitted 
into lead 3, we have here a — [3 = 7 = 3. 



t- 

CO 
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k B T = 0.1 

n = 0.5 

z = 0.3 
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slows, 
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r" C 3p. 



Vi 



FIG. 5: The third cumulant and its three contributions as a 
function of voltage. At high Vr the cubic term dominates. 
The dashed line Cg low shows the approximation in the high 
voltage regime. 



The radial direction can be interpreted as the energy of 
the particles. Due to the fluctuating potential, we allow 
energy exchange and average over all excitations during 
the traversal time r under the condition, that the energy 
takes the value E at point •. 

The modulation part Cf consists of three one-particle 
processes, that are connected via the correlation of the 
fluctuating potential only. 

In the second term C 3 S a modulated current and the 
noise are correlated also via the fluctuating potential. In 
the shot noise (the right ring in Fig. Efb)) we have two 
points X (time t) and • (time t + r), because we deal 
with a two-particle process: an electron propagates from 
channel p, to f3, a hole from (3 to v, an electron from v to 
7, and a hole from 7 to /i. In addition there is -even at 
zero temperature- a correlation between the modulated 
current and the Nyquist noise, which is not represented 
in Fig. El 

Cf p then is the actual three-particle exchange con- 
tribution, connecting the leads a, /3, and 7 by six S- 
matrices. Here we again average over excitations between 
the points •. 

In Fig. [5] the three contributions and the total third 
cumulant C 3 are depicted as a function of Vt. Note that 
C 3 is negative for $ = and z = 0.3, this implies that 
the distribution in lead 3 is skewed to the right. For 
high voltage the term C| p is linear, C 3 E is quadratic, 
and the cubic contribution Cf dominates. The full third 
cumulant is odd in voltage and thus vanishes for V = 0. 

In the symmetric case, 1Z = 0.5, we find AB- 
oscillations proportional to cos 3$ and cos$; with de- 
creasing z the terms in cos 3$ are stronger damped than 
those with cos$, as shown in Fig. Etc). Note also, that 
around $ = 0, ±7r the third cumulant changes sign with 
decreasing z. Fig. |3Id) shows the third cumulant as a 
function of temperature and voltage. In contrast to the 
noise the third cumulant contains no voltage indepen- 
dent (Nyquist) contribution whatsoever, and therefore, 
changes little with temperature^. 



We can distinguish two limiting regimes: high and low 
voltage^ in comparison to 1/r. In the first case, many 
particles (Vt 1) pick up approximately the same phase 
<p(t), since many wave packets fit at the same time into 
an interferometer arm. To the contrary for low voltage 
(Vt C 1) every wave packet acquires different phases 
tp(t), which are uncorrelated. Therefore, for high and 
low voltage we can make approximations of a slowly and 
a fast varying phase tp(t) respectively. 

We next compare these approximations with the exact 
expression. The modulation contribution Cf does not 
depend on the phase dynamics, but is determined by the 
parameter z. Therefore, it is independent of whether we 
consider a fast or a slowly varying phase. As Cf decreases 
with V 3 t 3 , in the low voltage regime it is small compared 
to the terms linear in voltage. In both limiting cases the 
averaging of the contributions C| p and Cj 3 s simplifies; ad- 
ditional details of the averaging and the approximations 
are discussed in the appendix B 

The low voltage approximation fits only for very small 
values of Vt, and we will not discuss this limiting case 
here further. A comparison of the high voltage approxi- 
mation, detailed in the appendix„ with the exact result is 
shown in Fig. [5] This approximation matches the exact 
curve very well for a wide range of Vt and even in the 
case Vt ~ 1. 



VI. GENERATING FUNCTION IN THE HIGH 
VOLTAGE REGIME 

The discussion above showed that in principle we can 
calculate any single cumulant with the only assumption 
of a Gaussian fluctuating potential. However the cal- 
culational effort grows very fast with the order of the 
cumulant. For this reason let us now concentrate on the 
limiting case of high voltage (Vt 3> 1), when many parti- 
cles pick up approximately the same phase. In this limit 
the generating function, Eq. (|SJ) can be evaluated directly. 
To this end, we perform the averaging, see Eq. JJjJ under 
the following two assumptions. 

(i) When the accumulated phase tp(t) = C +T dt'U(t') 
varies slowly, the matrix A(t, t') in Eq. JSJ takes the 
formic A(t, f ) -> n(t - t')(Si (<p)S(<p) - 1). It depends 
on a time difference and on a parameter ip, that itself 
is a function of a (single) time variable. After Fourier 
transformation we are left with one integral in time and 
another in energy. At zero temperature, the energy inte- 
gral is easily performed and one obtains 



Trln(l 



+A) ^i L° dtln[1 + T ^ melX - 1)] (14) 

Here Tsx(<p(t)) depends on the instantaneous phase ip(t). 

(ii) To take into account the correlation time r of 
phase fluctuations, we make an additional approxima- 
tion beyond the slow approximation discussed so far: We 
take the slowly varying phase <p(t) to be a steplike func- 
tion, that has constant values ip = IAt in time intervals 




FIG. 6: The distribution function (logarithmic scale) at 
zero temperature under influence of decoherence (N T = 10, 
to = t). (a) weak dephasing: z = e -0 ' 32 , (b) strong dephas- 
ing: z — e~ 2 . The flux dependence is suppressed, and the 
distribution tends towards a symmetric function for 1Z = 0.5. 



[fcr, (fc + l)r] for k £ N . The values ip and consequently 
the values U obey a Gaussian distribution. Note that the 
parameter U does not describe the actual fluctuations of 
the potential U(t) since these are determined by the cor- 
relation time trc which is here assumed to be smaller 
than t. 

The to jr time intervals of length r that can be treated 
independently. During each time interval the functional 
integral © becomes a simple Gaussian integral in the 
variable U : 



X(A) 



'dUe 



(l+T 31 (U)(e^ -1)) 



with T 31 (U) = 2TZT(1 + cos($ + Ut)). 
in TA leads to the generating function 



The integration 




JV T - 



™ a Pyrn-p z (2p-m) 2 



(16) 

with the constants c = 1 + 2RT(e lX - 1), 
a = KTe l *(e lX - 1), b = UTe"^ (e lX - 1). 

This approximation for the generating function applies 
also in the case trc > t: then the correlation time of the 
phase ip(t) is the same as the one of the potential U(t), 
trc, and the time r in Eq. 1|16[) has to be replaced by 
trc- 

Starting from this generating function (|16fl one finds 
immediately the cumulants by simply taking derivatives 
in A. These cumulants in their modulation contributions 
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from the slow limiting case of the cumulants presented in 
section|V] due to the treatment of cp(t) as a steplike func- 
tion (corresponding to a discretization in time space). 

Each cumulant shows harmonic oscillations in <E> up to 
the order of the derivative. As expected the terms con- 
taining oscillations with higher periodicity are damped 
stronger. Therefore, Eq. (|16|l is a convenient starting 
point to investigate the strong dephasing limit. In this 
case we can restrict the sum in Eq. (|lt>fl to small pow- 
ers in z. Up to first order in z the generating function 
consists of two hypergeometric functions F 



with 



F 



2 : 



Nr 

2 



1 



1 



ah 



2 ' '4c 2 



+c 



N T 1 N T 



ab 



(17) 
-.in. 



N T (c-l)zcos<£F(-— + -,-— + 1;2;-^ 



The first term represents the purely classical part, it does 
not depend on z and <E> but only on the beam splitter pa- 
rameter 1Z, A and N T . For TZ — 0.5 it is symmetric in A 
and thus yields a symmetric contribution to the distribu- 
tion function. The second part is the interference term 
proportional to z cos $; for TZ — 0.5 it is asymmetric in A. 

In the general case for arbitrary z, we perform the 
Fourier transform of Eq. i|15|) numerically to find the dis- 
tribution function P(Q). Fig. [5] depicts the logarithm of 
the distribution function for weak (Fig.UJa)) and strong 
dephasing (Fig. EOT))) for a time interval to = t. In the 
case of weak dephasing, it still exhibits the form of a bino- 
mial distribution with the periodicity in $ (compare with 
the no-dephasing case in Fig. Hl^a)). For increasing de- 
phasing strength, the AB-effect is suppressed. In the case 
of symmetric beam splitters the function P(Q) then tends 
towards a symmetric distribution around N T /2, and the 
third cumulant vanishes. However, for asymmetric beam 
splitters the distribution function is not symmetric. 

To extract information about the tails of the distribu- 
tion function P(Q), we use the stationary phase approx- 
imation and solve the saddle point equation d = iQ 
for A in the complex plane. For Q <C (Q), i.e. when 
nearly no particles are transfered to lead 3, we find 
(—A) ^> 1 for negative A. Starting from Eq. (|16fl we ob- 
tain 

lnf»(Q)«-Qln(— Q) + — lnu-Qln- (18) 

to T W 

The first term determines essentially the variation of the 
distribution function with Q. Only the parameter v and 
w are functions of the magnetic flux and the decoherence 
parameter z: 



rn—0p.q—0 

N T m 

w = J2 Vm pi 

rn—O p : q—0 



(19) 



N T -m + q 2TZF 
l~q 1 - 2KT 



N T \( rrt\f. N T —m\f m 



m/\p 



1 A? 



(20) 



Thus in Eq. (|18fl the dominating contribution of the Q- 
dependence of P(Q) is classical, and the last two terms 
represent the coherent part. For positive A 3> 1, the large 
Q tail, Q 3> (Q), of the distribution function is obtained 
in an analogous procedure. 



VII. CONCLUSION 

We analyze the Full Counting Statistics of a meso- 
scopic, electrical Mach-Zehnder interferometer pene- 
trated by an Aharonov-Bohm flux, and in the presence of 
a classical fluctuating potential. The Mach-Zehnder in- 
terferometer has the advantage that there are no closed 
orbits. Interference is only a consequence of superposi- 
tion of amplitudes in the out-going channels. Of interest 
is the generating function of transfered charge for such 
a quantum coherent conductor and the transition from 
coherent transport to classical transport. 

First, we discuss the probability distribution for a fully 
coherent Mach-Zehnder interferometer. At zero tempera- 
ture the distribution of transfered charge is binomial like 
in a two terminal conductor, but with a flux-dependent 
transmission probability. At finite temperature, the fact 
that we deal with a multi-terminal conductor leads to 
a different distribution which becomes classical due to 
flux-independent two-particle processes. 

In the presence of a fluctuating potential, assuming 
that it is Gaussian, we obtain the first three cumulants 
as a function of the Aharonov-Bohm flux and a dephasing 
parameter. The classical fluctuating potential modulates 
the conductance. As a consequence the noise contains a 
term proportional to voltage squared. The same effect 
results in a term cubic in voltage in the third cumulant. 
In the third cumulant the conductance modulation is cor- 
related as well with the shot noise. The cumulant of the 
k-th order shows oscillations in the flux up to cosfc$. As 
one expects the high-order harmonics of the Aharonov- 
Bohm oscillations are suppressed much stronger with in- 
creasing dephasing than the low order oscillations. 

For the limiting case of high voltage we find the gen- 
erating function. This permits, like in the coherent case, 
the calculation of the distribution function of the trans- 
fered charge. Therefore, in this limiting case, we are able 
to follow the suppression of the AB-oscillations in the 
distribution function. 

For the tails of the distribution function we find that 
they consist of two parts: a purely classical part, inde- 
pendent of dephasing, which determines the dependence 
on transfered charge and a quantum part which governs 
the flux dependence. 

The analysis of the conceptually simple case of a Mach- 
Zehnder interferometer will hopefully provide a useful 
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point of reference for future theoretical investigations of 
Counting Statistics in the presence of quantum fluctua- 
tions. However, the Mach-Zehnder interferometer is not 
only of theoretical interest, but recently has been realized 
experimentally 32 . The questions addressed here are thus 
also of interest for future experiments. 
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For the third cumulant four correlations functions 
h + ±± depending on two time variables are possible, for 
example 



fc + _+(t,f) = i ((e 



z)(e 



-iip(t') 



z)(e 



Mo) 



-z) 
(A.3) 

The functions h(t, t') can be expressed in terms of 
g± given by Eq. (|A.2|1 . 

With Eqs. and 1A.1|) all cumulants take a compact 
form. The mean current in lead 3 is 



= — Tr[n(K) - nP], (A.4) 

A=0/ U *° 



APPENDIX: CALCULATION OF SINGLE 
CUMULANTS 



leading to the result presented in section Ivl 

For the noise we express the second cumulant in terms 
of derivatives of A, taken at A = 0: 



Successive derivatives in A of the logarithm of the gen- 
erating function, see Eq. yield the cumulants. Using 
the expression JSJ each cumulant is expressed in terms of 
derivatives of the matrix A, introduced in section I1VI 
Derivatives of A of A:-th order can be written in terms 
of a matrix K = PS and the matrix n, defined in sec- 
tion llVI P is a projection matrix with only non- vanishing 
element, P 33 = 1. Note that K has non- vanishing ele- 
ments only in the upper diagonal 2x2 block. Because 
backscattering is excluded in the MZI, the product PK 
vanishes and we obtain the simple form 



d k A 



dX 



i k n 



(K+(-l) k P) 



(A.l) 



The matrix K is separated into constant matrices (K) 
and L, and time-dependent scalar functions as follows: 
K = (K) + L(e^W - z) + Z,t( e -i<pW - z ). Here z is the 
decoherence parameter z — ^e ±lip< - t) ) u , introduced in sec- 
tion The average is calculated by replacing the gen- 
erating function in Eq. (JSJl by the phase factor; at this 
point we explicitly assume a Gaussian fluctuating poten- 
tial U(t), and use trc "C t. 

The average of the matrix K is denoted by (K), it 
is determined only by the parameters 1Z, and z. 
The averaging of the cumulants only affects the time- 
dependent functions (e ±lv ^ — z): Depending on the or- 
der k of the cumulants one obtains products of the func- 
tions (e ±lv W — z) at k different times. The averaging 
then defines correlation functions of different order. 

In the noise we get two correlation functions, depend- 
ing on one time variable and determined by the parame- 
ter z: 



. r„±2(l-M) _ 



- •>>, 



l]e(r-|*|). 



(A.2) 



The procedure of the average is similar to the calculation 
of the parameter z. 



JL fr _/t — \ 
iH Q \ \ 1 dX \ r dX J 



1 d 2 A 
+ W \ dA 2 " 



Tr 



dA 
~dX 




(A.5) 



The first term corresponds to the modulation noise Si 
and disappears in the case of no dephasing. The second 
term is the two-particle contribution (that separates into 
S nyq and S ox as explained in section^. Introducing the 
if-matrix, Eq. (|A.5|) simplifies to 

E 



U(jT[nK-n{K)]f 

+ — ( r Tr[nK(l - n)K] 
to \ 



(A.6) 



since the product K (1 — K) vanishes. The final result 
contains the Fourier transforms of the functions g± and 
a frequency integral, as also found in Ref. [33]. 

As for the current and the noise, we express the third 
cumulant in terms of derivatives of A at A = 0: 



dA I dA 
Tl dX \ Tl dx 





(A.7) 



d?A (dA 

Tr dx^- Tl (dX 



, dA d 2 A 
~ 3 ' l : ^ d* 




The first term corresponds to the modulation contribu- 
tion Cj 3 , the second term to the term C 3 E , and the re- 
maining part gives the three-particle contribution C| p . 
One checks immediately, that the first two terms vanish 
in the absence of dephasing. Using again the compact 
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expressions in K, Eq. (|A.7(1 becomes 



C 3 = - 



(Tr[nK-n(K)}) 3 ^+ (A.; 
+3 l(Tr[nK - n (K)]) [Tr[nK(l - ri)K})\ + 
+ ( r Tt[nK{l-n)K{l-2n)K 1 \ 



The evaluation of Eq. (|A.8() is quite lengthy. Introducing 
Fourier transforms of the correlation functions h and g 
and of the occupation functions n describing the leads, 
one is left with an integral over two frequencies. This 
integral contains a convolution of h or g with temper- 
ature and voltage dependent kernels. We evaluate the 
integration numerically. 

The limiting cases were introduced in section[V] In the 
high voltage regime, the accumulated phase varies slowly. 
Therefore, we evaluate the correlation functions g± and 
h (see Eqs. (JA.2(1 and (|A.3(l ) at equal times: 



slow: 



9±(t) 
h(t,t') 



fl±(0) 
h(0,0). 



(A.9) 



In this limit the correlation functions are entirely deter- 
mined by the parameter z. The frequency integrals in 



the final expression of C| p are given by integrands at 
zero frequency and thus simplify substantially. As the 
modulation term Cf does not depend on the phase dy- 
namics, it does not change in this approximation. The 
term C 3 S combines a modulation part and a two-particle 
process. Here we use h(t,t') w h(0,t'). 

The high voltage approximation does not catch contri- 
butions to the third cumulant due to particle-hole exci- 
tations created by the fluctuating potential U{t). Simi- 
lar to Ref. [33], this leads to an unimportant linear off- 
set between the exact result and the high voltage result. 
For better comparison, this offset is subtracted in Fig. EI 
which shows a good agreement of the high voltage ap- 
proximation with the exact result. 

In the opposite case, for low voltage, the phase factors 
in formula (|A.2|) and (|A.3(1 at different times are uncorre- 
cted and thus can be averaged independently and vanish. 

-t3 
y 3p 



In this case, we use in Cf„ the approximation 



fast: 



9±(t) 
h(t,t') 



(A.10) 



C 3 is proportional to V 3 r 3 and thus in the low voltage 
regime small compared to the terms linear in voltage. For 
low temperature also the contribution Cj 3 s is negligible. 
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